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We investigate the magnetic instabilities of the two-dimensional model of interacting e^, electrons for hole 
doping away from two electrons per site in the mean-field approximation. In particular, we address the oc- 
currence of orbitally polarized states due to the inequivalent orbitals, and their interplay with ferromagnetic 
and antiferromagnetic spin order. The role played by the Hund's exchange coupling Jh and by the crystal 
field orbital splitting Ez in stabilizing one of the competing phases is discussed in detail. 
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1 Introduction The richness of cooperative phenomena encountered in transition metal oxides con- 
tinues attracting considerable attention. In addition to the spin and charge degrees of freedom, also the 
coupling to the lattice and, in particular, the orbital degrees of freedom at and close to the orbital degen- 
eracy lead to very interesting behavior Here we will focus on the phenomena observed in nickelates 
and manganites, which can be attributed to strongly correlated Cg electrons. So far, it is known that in the 
regime of large intraorbital Coulomb interaction U strong quantum fluctuations may lead to qualitatively 
new behavior in a Mott insulator with three Cg electrons per site |2|, but the competition between various 
magnetic and orbital instabilities was little explored in the weak coupling regime. 

While several features are generic in the models with two f3l or more 14] orbitals per ion, and occur 
already when diagonal hopping is assumed, we note that these models are closer to the behavior of t2g 
electrons — strong interactions between them might explain a ferromagnetic (FM) instability in ruthenates 
f5]. In contrast, the orbital flavor for eg electrons is not conserved, and this is likely to lead to partial 
orbital polarization which is expected to modify the magnetic instabilities. Here we will consider such a 
realistic two-dimensional (2D) model of eg electrons with (ddcr) hopping element t at intermediate and 
strong coupling, which includes the full structure of on-site Coulomb interactions with two parameters: 
the Hubbard element U and Hund's exchange Jh |6|, and the crystal field orbital splitting E^- 

The magnetic and orbital instabilities within the eg band were less investigated until now, but they 
become very relevant in the context of doped La2-a;Srj:Ni04 nickelates, where interesting novel phases 
including the stripe order were discovered |7). They were obtained in the theory using realistic models 
both in Hartree-Fock ijSj and in exact diagonaUzation of finite clusters including the coupling to the lattice 
f9]. Furthermore, the role of Hund's exchange in the FM instability and in the metal-insulator transition 
was emphasized using a multiband Gutzwiller wave-function |10|. All these studies reveal interesting 
competition between FM and antifeiTomagnetic (AF) instabilities at density n = 1, and the competition 
between the latter two and a C-AF instability at n = 1.5. In this paper we investigate the pure electronic 
problem for eg electrons, using the mean-field approximation, seeking for phases which are both orbitally 
and magnetically polarized. By varying the electron density n between half filling n = 2 and n = 0, we 
cover the hole doping regime a: = n — 2 relevant to the nickelates. 
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2 Model In this work we consider a two-dimensional (2D) model for eg electrons on a square lattice, 

H = Hkin + Hint + Hz, (1) 

with two orbital flavors: \x) ^ \x'^ — y^) and \z) ~ |3z^ — r^). The kinetic energy is described by 

{ij) al3a ^ ' 

where t stands for an effective [dda) hopping matrix element due to the hybridization with oxygen orbitals 
on Ni— O— Ni bonds, and the off-diagonal hopping ifj along a and h axis depends on the phase of the |a;) 
orbital along the considered cubic direction. The electron-electron interactions are described by the on-site 
terms, which we write in the following form l6l . 

i i 
i i 

with rijQ, = ^Ylia n-iaa (for a = X, z). U and Jh stand for the intraorbital Coulomb and Hund's exchange 
elements. The interactions Hint are rotationally invariant both in the spin and in the orbital space. The last 
term H^ describes the uniform crystal-field splitting between — y"^ and — orbitals, 

Hz^ ^EzY,i'i,,-n^z). (4) 

i 

It is convenient to rewrite the Hamiltonian (|3} by introducing the following operators: 

rii = ^^fliaa, (5) 

aa 

TTli — ^ ^ ^(T^iaa : ^ ^ ^a^iacr : Pi — ^ ^ ^a^a^^iaa : (6) 



fia — (^iaa'^a0'^i0a' 0) 

a/3 

with Act — ±1 for a =t (i) spin and Aq = ±1 for a — x{z) orbital, and cr^ is a Pauli matrix. These 
operators correspond to the total density, the total magnetization, the orbital polarization, the magnetic 
orbital polarization, and the on-site orbital flip respectively. The Coulomb interaction term Hint (E) can be 
then written as: 

i 

+ JHY^hhi- (8) 

i 

The order parameters introduced in Eqs. (|6|l: m^, Oi, and pi, used next to minimize the ground state 
energy, provide a complete description of the ground state at finite doping. We emphasize that they also 
reveal the dominating role of the kinetic energy of doped holes over the superexchange energy. Namely, 
large electron filling of | z) orbitals, contributing to a naiTow band, optimizes the kinetic energy of holes in 
magnetically polarized states, doped into the \x) orbitals, contributing to a wide band. On the contrary, the 
superexchange oc J = A{t"")'^/U at large U suggests that the system would better optimize the magnetic 
energy when the orbitals with larger hopping elements t^^ are closer to half filling. We show below that 
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Fig. 1 Order parameters: magnetization (m), orbital polarization (o), and magnetic polarization (p) in the FM 
phase as a function of doping x for: (a)-(c) Jh = 0.25f/ and (d)-(f) Jh ~ 0.1517, and for two values of the Stoner 
parameter: U + Jh ~ 8t (black solid lines) and U + Jh = 4f (gray solid lines). A and B refer to two sublattices in 
the orbital ordered state for x ^ 1. The orbital polarization in the reference PM states is shown by dashed lines. 



the complex interplay between all the degrees of freedom of the model Q results in rather peculiar doping 
dependence of the order parameters (|6jl, and thus leads to highly nontrivial and rich phase diagrams. 

We investigated the stability of possible phases with either uniform or staggered magnetic order in the 
mean-field (Hartree) approximation by expressing the local operators (|6j by their mean-field averages, 
jf ~ 2ji{'ji) — (7i)^. In order to establish unbiased results, the calculations were carried out on a large 
128x128 cluster, using periodic boundary conditions at low temperature T — Q.Olt (here fcs — 1). 
Consistently with the present mean-field analysis we assumed {fi„) = 0. 



3 Numerical results 



3.1 Magnetic order and orbital polarization: First of all, in the paramagnetic (PM) state at = 0, 
with (m) = (p) = 0, a higher electron density is found in \x) orbitals ((o) > 0), as then the kinetic energy 
is lowered [Figs.^b) anddJe)], except at a; = 0. This state is our reference state for possible magnetic 
instabilities. 

We now proceed with the discussion of the magnetic order and orbital polarization for two characteristic 
values of the Stoner parameter I = U + Jh- intermediate coupling / — At, and strong coupling / — 8t, 
being smaller and larger than the bandwidth W = 6t, respectively. Let us begin with the FM phase. In 
Fig.^a) we show the magnetization (m) as a function of doping x for the ratio Jh /U = 0.25 which is 
representative of the strong Hund's exchange coupling regime. In this case the interaction in the o-channel 
is repulsive. As shown in Fig. [J several FM phases occur. 

Consider first the intermediate interaction strength I = At [Figs.^a-c)], where one finds two discon- 
nected FM states: one for < a; < 1, and the second one for x ~ 1.5. The latter corresponds to a van 
Hove singularity in the density of states. Since it is predominantly related to the \x) orbital, both the orbital 
polarization (o) [Fig.Qlb)] and the magnetic polarization (p) [Fig.Qlc)] are positive in this doping regime. 
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Fig. 2 Order parameters as in Fig.Qbut for the AF phase. 



In contrast, for < a; < 1 the total energy is minimized when a higher electron density is found in the 
\z) orbital, where also larger magnetic moments are formed. Such anisotropic filling of eg orbitals follows 
from a large difference between the t^^ and t^^ hopping elements 1 1 1 1. Note that the orbital polarization 
is here opposite to that in the reference PM state. 

The above peculiar behavior disappears gradually when the interaction strength is enhanced to / = 8i 
and both (o) and (p) tend to saturate to the optimal value, being positive (negative) for a; < 1 (a; > 1), 
leaving the \z) orbital almost fully polarized. Here the magnetic instability, with the largest effective 
interaction term \J + Ju, dominates and the magnetization barely deviates from its saturation value 2 — 
except for the low electron density a; > 1.8. In this case doping the half-filled FM state first leads to holes 
introduced into the \x) orbitals, leaving the center of the narrower band, with predominantly \z) orbital 
character, below the former broader one. Therefore, in this situation the formation of localized magnetic 
moments optimizes the energy. They are naturally associated with the \z) orbitals since they contribute to 
the naiTower band. 

When Hund's exchange coupling Jh is reduced, the interaction in the o-channel becomes attractive. As 
a result, for large / = Si, both (o) and (p) nearly saturate to the ideal behaviors —a; for a; < 1 and 2 — x for 
a; > 1. Therefore, a transition between these two solutions would be first order, and one observes a jump at 
a; ~ 1. However, as shown in Figs.[3d-f), the two-sublattice orbital order sets in in this crossover regime in 
the form of a FMj;^ state. This state has opposite orbital polarization (o) a ~ —0.8 and (o) b ~ 0.8 on both 
sublattices [Fig.ffle)]. While the total density (ri) a is somewhat higher than (rt) b due to inequivalent Cg 
orbitals, also (m) a > {m)B- For large doping a; > 1 the electrons occupy mainly the |a;) orbitals, and the 
small occupancy of \z) orbital results solely from the interorbital hopping term cx t^z- Indeed, at a; 1.3 
one finds appreciable orbital polarization, with |a;) orbitals occupied and almost empty \z) orbitals, the 
situation encountered in Lai_a;Sri+2;Mn04 manganites 1 12|. In all FM phases found at Jh = 0-15U the 
total magnetization is close to saturation. When U is reduced, one gradually recovers the behavior obtained 
for large Jh/U. 

We now turn to the AF phase, expected as a ground state near half filling (a; — 0). The order parameters 
are shown in Fig.|2]for the same parameter values as for the FM case. For I = 8t and Jh/U — 0.25 the 
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Fig. 3 Order parameters as in Fig.Qbut for the C-AF phase. 

mean-field equations possess two competing solutions [Figs.|3a-c)]. The first one, which can be continued 
to weak coupling, is characterized by negative values of orbital (o) and magnetic (p) polarizations. Namely, 
the higher electron density is found within the \z) orbitals, and these orbitals carry the magnetic moment. 
More precisely, introducing holes in the half-filled insulating AF state mostly affects the band with |a;) 
orbital character, leaving the localized magnetic moments within the \z) orbital almost saturated. This 
solution extends to large doping x ~ 1. In contrast, the second solution rather stems from the behavior 
expected for low density: the electrons first occupy the broader band with |.t) orbital character until quarter 
filling {x = 1) is reached, and next they gradually occupy the other band, with \z) orbital character. 
However, since the interaction in the o-channel is repulsive and since both bands are coupled, the orbital 
polarization (o) and the magnetic polarization {p) are reduced from their maximal values which would 
be reached for decoupled orbitals. These two trends from a; = and x — 2 contradict each other, and 
therefore an abrupt (first order) transition between both solutions is observed at a; ~ 0.7 [Figs.|2jb-c)]. 

Reducing Jh barely affects the above findings for strong coupling I — 8t [Fig. |2jd-f)]. While the 
magnetization (m) is almost unchanged, the first order transition between two differently polarized states 
is more pronounced, as the values of the orbital polarization and the magnetic polarization are enhanced. 
When reducing U the location of the first order phase transition shifts towards smaller doping, while all 
order parameters are suppressed. At the same time the critical doping locating the second order phase tran- 
sition is reduced by a weaker Coulomb interaction U but is enhanced by a weaker Hund's exchange cou- 
pling Jh- When seeking for other phases one may expect that two-sublattice FM solutions can smoothly 
interpolate between the FM and AF states. Such solutions never turned out to be the ground state in this 
study (up to / = 8t). 

A competition between the FM and AF order in the present model of Cg band may lead to a superposition 
of the two phases in a form of C-AF phase, where the magnetic moments are FM along one direction and 
staggered in the other (orthogonal) one. According to recent numerical simulations |9|, a coexistence of 
FM and AF bonds is indeed expected for x ~ 0.5. Unlike in the AF phase, the order parameters are 
continuous functions of doping for Jh /U = 0.25, as can be seen in Figs.|3la-c). This behavior is similar 
to that of the FM case (Fig.Q. Its origin can be attributed to the orbital polarization which is substantially 
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Fig. 4 Phase diagrams of the Cg orbital model Q as functions of the Stoner parameter U + Jh and hole doping 
X = 2 — n, with: (a) Jh ~ 0.2517, and (b) Jh = 0.15t/. Panel (c) shows the stable phases for finite crystal field 
splitting Ez = 2t and Jh ~ 0.25U. Transitions from the PM phase to magnetic phases are second order. The 
remaining solid lines denote first order transitions while the dashed, dotted, and dashed-dotted lines indicate second 
order transitions. 



stronger in the AF case to the extent that is exceeds a certain threshold above which no smooth solution 
can interpolate between the small and large doping regimes. For x ~ 0.5 the magnetic moment (m) is 
caiTied by the \z) orbital for large U, while (m) decreases and (o) changes sign for small U. 

When reducing Jh /U, the orbital polarization is enhanced and a first order transition appears for I = 8t 
[see Figs. |3ld-f)]. In this case both the orbital polarization and the magnetic moment are predominantly 
carried by the stronger correlated \z) orbital (with a weaker hopping and thus larger ratio U/t^^ than 
U/t^^) in the physically relevant doping range centered around x = 0.5. 

3.2 Magnetic phase diagrams: Our main findings are summarized in the phase diagrams in Fig. |4] For 
Jh /U = 0.25 [Fig.01a)], the doped PM phase is characterized by a positive orbital polarization, therefore 
denoted PMx. It is unstable towards AFx phase for small doping up to x ~ 0.5, towards C-AFx phase 
for 1 < X < 1.02 and for 1.65 < x < 1.75, and towards FMx phase otherwise. In particular, for 
X = 1 the FMx, C-AFx, orbitally unpolarized FM and (for / > 16Q the alternating FMxz phases appear 
successively with increasing interaction strength /. The C-AF phases are found for x ~ 0.5 at / > At, and 
also in a small range around x ~ 7/4. 

When reducing Jh/U [Fig. |4](b)], the main difference appears for x ^ 1. Here the PMx phase is 
unstable towards an AFx phase which itself is robust and remains stable up to strong coupling. This 
seemingly peculiar behavior can be better understood by diagonalizing exactly a two-site cluster. One can 
find the ground state to be AFx for small Jh and arbitrary U and FMxz for large Jh/U, in qualitative 
agreement with our mean-field calculation. 

Let us finally mention that the model we use is known to have a van Hove singularity in the vicinity 
of X = 1.5, which is expected to induce a FM instability for arbitrary weak coupling at zero temperature. 
This particular instability, however, turns out to be unusually strongly temperature dependent. Therefore, 
the corresponding critical value of the Stoner parameter / is finite at temperature T > 0, and reaches a 
value close to W/i (Fig. at T = 11^/600 used in this work. 

3.3 Consequences of the crystal field splitting: A complete investigation of the phase diagrams at finite 
crystal field splitting Ez would be quite involved, and is left for future work. Ai Ez = the majority of 
stable magnetic solutions is characterized by a positive orbital polarization, a tendency expected for a 2D 
model of tg electrons 1 12] , which would certainly be enhanced by a negative crystal field Ez- We therefore 
limit our present discussion to the influence of a positive Ez, in order to investigate a competition between 
the kinetic energy, which is lower when the broad band with predominantly |x) orbital character is closer 



© 2003 WILEY- VCH Veikg GmbH & Co. KGaA, Wcinhcim 



pss header will be provided by the publisher 



7 



to half filling, and the potential energy at finite Ez- Furthermore, we restrict ourselves to the strong Hund's 
exchange coupling regime Jh ~ 0.25[/. 

As shown in Fig. |3c), the orbital polarization of the PM phase is changed to negative (PMz phase) 
already for moderate = 2t. As a result, the magnetic moment and the orbital polarization are carried 
by the same orbital in all phases, and the magnetic instabilities are achieved for lower values of / = 
U + Jh- Another consequence of finite E^ > is the observed shift of the van Hove singularity to 
larger doping, strongly enhancing the tendency towards ferromagnetism in the low density regime. In 
addition, the competition between FM and AF phases at quarter filling (x = 1) remains quite spectacular: 
even though both weak coupling and strong coupling expansions predict antiferromagnetism, as does our 
calculation, FM phase nevertheless takes over in an intermediate coupling regime 5.15t < / < 6.5i. 

Such features as the AFz phase obtained for low doping, the C-AFz phase for x ^ 0.5, and FM one for 
both ~ 0.5 < a; < 0.98 and x > 1.02, are robust and are found also at finite crystal field, in particular in 
the strong coupling regime. On the contrary, the instability of the PM phase to the FMa; phase disappears, 
and the one to the C-AF phase moves from x ^ 1.5 to small doping. 

4 Conclusions In summary, we have determined the phase diagram of Cg electrons on the square lattice 
within the mean-field approximation. The occurrence of antiferromagnetism in the vicinity of half filling, 
followed by C-AF phase at a; ~ 0.5, as well as FM phases for x ~ 0.75 and a; ~ 1.5, are robust features of 
this model. Note that the regions of stability of the AF and C-AF phases with respect to the FM one would 
still be somewhat extended due to quantum corrections 1131 . In particular, the occurrence of C-AF phase 
indicates that even more complex types of magnetic phases, such as stripe phases with larger magnetic unit 
cells 1 7 9 1, might be expected in doped nickelates. 

In contrast, the ground state is strongly parameter dependent in the vicinity of quarter filling, x — 1, 
resembling to some extent the ground state of the model with two equivalent orbitals |3 5|. While the 
orbital polarization systematically appears in all phases, the orbital carrying the magnetic moment does 
not necessarily coincide with the one carrying higher electron density, leading to a particularly interesting 
interplay between magnetic and orbital degrees of freedom. 
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